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We present the discovery of OB type absorption lines superimposed to the emission 
line spectrum and the first double-lined orbital elements for the massive Wolf-Rayet 
binary HDE 318016 (=WR 98), a spectroscopic binary in a circular orbit with a 
period of 47.825 days. The semiamplitudes of the orbital motion of the emission lines 
differ from line to line, indicating mass ratios between 1 and 1.7 for Mwr/Mob- 
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1 Introduction 



HDE 318016 was discovered to have an emission line spectrum of Wolf-Rayet 
type by Cannon Sz Mayall (1938). Because both N and C emission lines appear 
strong in the spectrum it was classified as WC7-N6 by Smith (1968). This star 
was included in the Sixth Catalogue of Galactic Wolf-Rayet Stars (van der 
Hucht et al. 1981) as WR 98 and classified as WN7+WC7. However WR 98 
was confirmed to be a single-lined binary with a period of 47.8 days with N 
and C emission lines moving in phase (Niemela 1991). 

The optical spectrum of WR 98 has been described by Lundstrom & Stenholm 
(1984). Conti & Massey (1989) proposed the nomeclature WN/WC for WR 98 
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suggesting that this star, along with others with similar emission line spectra 
where both N and C lines are observed, has a transition composition between 
the WN and WC subclasses. 

Smith et al. (1996) classified WR 98 as WN8o/C7, using a new three-dimensio- 
nal classification scheme, where the "o" means that no hydrogen is observed 
in the spectrum. Relevant parameters of WR 98 can be found in the recent 
Vllth Catalogue of galactic WR stars (van der Hucht 2001). 

WR 98 is a probable member of the open cluster Trumpler 27 (e.g. Feinstein 
et al. 2000), and it has also been detected as a non-thermal radio source by 
Abbott et al. (1986). Furthermore, WR 98 exhibits random relatively large 
amplitude optical light variations (~ 0.1 magnitude in the Johnson V filter), 
typical of stars with WN8 type spectrum, but a periodicity for these variations 
has not been found (Marchenko et al. 1998). 

In this work we present a detailed radial velocity analysis of optical spectral 
lines of WR 98 showing it to be a double-lined binary with high minimum 

masses. 



2 Observations 

2.1 Photographic spectra 

41 photographic spectrograms were obtained at the Cerro Tololo Inter- Ameri- 
can Observatory, Chile, between 1980 and 1984. These spectrograms were 
secured with the image-tube spectrograph at the Cassegrain focus of the 1-m 
Yale telescope. All exposures were made on Kodak Illa-J emulsion baked in 
"forming gas" (N 2 + H 2 ). The spectrograms have a reciprocal dispersion of 45 
A mm" 1 . A spectral region from ~ 3600 to 5000 A was covered. Exposure 
times vary between thirty and ninety minutes, giving Signal to Noise (S/N) 
ratios ~ 15-30. A He-Ar lamp spectrum was used as wavelength calibration 
source. 

A preliminary report of the binary nature of WR 98 based on these data was 
presented in the IAU Symp. 143 (Niemela 1991). 

2.2 Digital spectra 

28 optical digital spectral images of WR 98 were obtained with the Cassegrain 
Boiler & Chivens (B&C) and REOSC spectrographs attached to the 2.15-m 
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telescope at Complejo Astronomico El Leoncito (CASLEO) 3 in San Juan, 
Argentina, between 1997 and 2001. Nine spectra were secured with the B&C 
spectrograph, in 1997, March, and 1998, February and May. A PM 512x512 
pixels CCD, with pixel size of 20/zm, was used as detector. The reciprocal 
dispersion was ~ 2.3 A pixel -1 , and the wavelength region covered was about 
AA 3900 - 4900 A. 

Nineteen spectra were obtained with the REOSC spectrograph, between 1999, 
March, and 2001, October. For these spectra a TEK 1024x1024 pixels CCD, 
with pixel size of 24/im, was used as detector. The reciprocal dispersion was 
~ 1.6 A pixel -1 , and the observed wavelength region was ~ AA 3850 - 5450 A. 

We used a slit width of 2 arcsecs for all our spectra. Typical exposure times 
for the stellar images were between 30 and 40 minutes, resulting in spectra of 
signal-to-noise ratio S/N ~ 50-100. 

He-Ar (or Cu-Ar with REOSC spectrograph) comparison arc images were 
observed at the same telescope position as the stellar images inmediately after 
or before the stellar exposures. Also bias and flat-field frames were obtained 
each night, as well as flux and radial velocity standard stars. 

All spectra were processed with IRAF 4 routines at La Plata Observatory. 
2. 3 Determination of radial velocities 

Radial velocities of lines in the photografic spectra of WR 98 were measured 
with the Grant oscilloscope comparator-microphotometer at the Instituto de 
Astronomia y Fisica del Espacio (IAFE), Buenos Aires, Argentina. 

For the present study, selected photografic spectrograms and their correspond- 
ing calibrations were digitized with a Grant microphotometer, and calibrated 
and measured with IRAF routines. 

The emission lines in the spectrum of WR 98 appear approximately gaussian 
in shape (cf. Fig 1). Therefore, for determination of radial velocities in our 
digital spectra, we measured central wavelengths of all emission lines fitting 
gaussian functions to the line profiles. 

The journal of the spectroscopic observations with the radial velocity measure- 
ments is presented in Table 1, where succesive columns quote: the Heliocentric 

3 operated under agreement between CONICET, SeCyT, and the Universities of 
La Plata, Cordoba and San Juan, Argentina 

4 IRAF is distributed by the NOAO, operated by the AURA, Inc, under cooperative 
agreement with the NSF, USA. 
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Julian Date of each observation, the orbital phase, and the radial velocities 
for Niv, Nv, Cm, Hell emission lines and H<5 and H7 absorption lines, respec- 
tively. 



Table 1: Observed heliocentric radial velocities (in km s x ) in the spectra of 
WR 98 



HJD 


Phase 1 


N IV em. 


N V em. 


Cm em. 


Hen em. 


H7 abs. 


H & abs. 
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44744.820 


0.521 






-79 


-38 


43 




44895.510 


0.672 


-118 


-57 


-93 


-74 






44897.520 


0.714 


-188 


-171 


-102 


-44 


160 




44898.520 


0.735 


-226 


-60 


-113 


-106 






45067.880 


0.276 


-4 




11 


38 


-47 




45068.890 


0.297 


-5 


87 


33 


42 


-128 




45069.910 


0.319 




79 


46 


88 






45070.910 


0.339 


-27 


73 


-9 


20 


-103 




45123.840 


0.446 


88 




-10 


30 






45184.640 


0.718 


-183 




-110 


-47 






45186.590 


0.758 


-190 




-142 


-88 






45189.610 


0.821 






-96 


-63 






45190.600 


0.842 




-71 


-125 


-67 






45191.620 


0.863 


-65 


-166 


-65 


-49 






45192.630 


0.885 


-208 


-109 


-79 


-16 






45196.630 


0.968 


-51 


-76 


-77 


-24 






45197.670 


0.990 


-59 


57 


-78 


-19 






45198.640 


0.010 


-79 


80 


-53 


11 






45199.630 


0.031 


-29 


29 


-16 


9 






45248.520 


0.053 


-25 


55 


-35 


25 




-221: 


45249.520 


0.074 


36 


121 


21 


43 


-46 




45250.510 


0.095 


83 


64 


76 


86 


-115 




45251.540 


0.116 


93 


80 


-1 


68 




12 


45252.540 


0.137 


9 


37 


40 


66 


-80 


-27 


45254.530 


0.179 


84 


131 


-17 


30 






45255.530 


0.200 


58 


-17 


-12 


23 




-72 


45256.530 


0.221 


139 


96 


76 


113 
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Table 1 Continued 



HJD 

2 400 000+ 


Phase 1 


N IV em. 
A 4057.76 


N V em. 
A 4603.73 


Cm em. 
A 4648.83 


Hen em. 
A 4685.68 


H7 abs. 
A 4340.47 


H (3 abs. 
A 4861.33 


45257.540 


0.242 




-26 


74 


102 






45258.530 


0.263 




16 


21 


48 






45260.530 


0.304 


-5 


99 


50 


29 




-67 


45502.750 


0.369 




113 


-34 


36 






45506.840 


0.455 


68 


-13 


5 


62 






45507.830 


0.475 





37 


-19 


40 






45508.780 


0.495 




-31 


-11 


24 


6 


92 


45555.670 


0.476 


10 


-57 


-39 


4 






45846.750 


0.562 


12 


7 


-61 


-37 




206: 


50536.880 


0.631 


-120 


-51 


-103 


-41 


116 


88 


50537.883 


0.652 


-123 


-133 


-109 


-53 






50538.872 


0.672 


-124 


-118 


-112 


-34 




118 


50539.886 


0.693 


-125 


-104 


-122 


-78 




39 


50540.878 


0.714 


-176 


-144 


-130 


-55 




86 


50541.882 


0.735 


-130 


-143 


-136 


-59 






50542.877 


0.756 


-169 


-141 


-100 


-70 




144 


50858.851 


0.363 


-38 


50 


-34 


26 


-24 


-120 


50860.844 


0.404 


-45 


5 


-47 


15 


-47 


-124 


51261.900 


0.790 


-93 


-61 


-121 


-41 






51261.913 


0.791 


-125 




-125 


-38 






51264.909 


0.853 


-106 


-75 


-122 


-49 






51266.907 


0.895 


-96 


-4 


-104 


-31 






51267.901 


0.916 


-37 


-28 


-91 


-34 




39 


51303.768 


0.666 


-102 


-126 


-105 


-47 


50 


135 


51304.820 


0.688 


-192 




-112 


-41 






51305.814 


0.709 


-160 


-136 


-119 


-46 




184 


51653.819 


0.985 


-86 


-17 


-68 


17 




-15 


51654.833 


0.006 


-69 


-37 


-56 


12 




47 


51806.573 


0.179 


134 


110 


8 


110 




-191 


51807.534 


0.199 


87 


66 


29 


81 


-68 


-189 


52009.863 


0.430 


-5 


-14 


-27 


47 


-57 


-74 


52010.898 


0.452 


-21 


-22 


-2 


42 




-7 
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Table 1 Continued 



HJD 

2 400 000+ 


Phase 1 


N IV em. 
A 4057.76 


N V em. 
A 4603.73 


Cm em. 
A 4648.83 


Hen em. 
A 4685.68 


H7 abs. 
A 4340.47 


H (3 abs. 
A 4861.33 


52011.880 


0.472 


-15 


-15 


-41 


37 


42 




52192.503 


0.249 


51 


36 


18 


52 




-169 


52193.496 


0.270 


69 


90 


29 


83 




-68 


52194.503 


0.291 


66 


69 


16 


64 




-97 


52195.501 


0.312 


52 


73 


7 


59 




-161 



1 Phases were computed according to To = 2,445,676.4+47.825E (See below). 



3 Results and Discussion 

3.1 The spectrum 

The blue optical spectrum of WR 98 is illustrated in fig 1, with identifica- 
tions for main spectral features. As seen in fig 1, the spectrum of WR 98 is 
dominated by emission lines of Nlll, Niv, Hell, Hel, with a strong Cm feature 
at A 4650A , consistent with the WN-WC classification originally proposed by 
Smith (1968). All of Hel emission lines in our blue spectra show P-Cygni pro- 
files, as well as Hell 541lA and Niv 5203A (See Fig. 1). Relative intensities 
of Nlll, Hell, Niv, and Cm lines in our spectra of WR 98 indicate a spectral 
type WN7-8/C. No hydrogen is detected in the spectrum. Relative intensities 
of emission lines in the spectrum of WR 98 show no appreciable changes in 
the whole of our dataset spanning 20 years of observations. 

We have compared WN the spectrum of WR 98 with that of stars classified as 
WN7o and WN8o (Smith et al. 1996), namely WR 55 and WR 123, for which 
we also have digital spectra observed with the same instrumental configuration 
as WR 98. This comparison indicates that WR 98 has a higher ionization 
degree than the WN8o star, because Niv emission lines in the spectrum of 
WR 98 are stronger and Nm emission lines are weaker. Actually, the WN 
spectrum of WR 98 mostly resembles that of WR 55, classified as WN7 by 
Smith et al. (1996). Thus our blue optical spectra of WR 98 would be best 
described as WN7o/WC. Fig. 2 shows the spectra of WR123, WR98, and 
WR55 for comparison. 

Several faint absorption lines were detected upon the WN emission lines in 
our spectra of WR 98. These absorptions are H7 and H/3, Hel AA 4026, 4471, 
and 5015 A and Hell AA 4200, and 4540A . As will be shown below, these 
lines belong to an OB companion. Because Hell absorptions appear fainter 
than Hel in the OB spectrum, we presume that the companion probably is of 
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4000 ^~ ~~ ' 4400 ^ ~~ ' 4800 

Wavelength [1] 



Fig. 1. Continuum normalized spectrum of WR98 obtained at CASLEO in 2000, 
September. Main emission lines are identified above, and companion OB absorptions 
below the spectrum 

spectral type 08-9. The luminosity class is not possible to determine from our 
data. 



3.2 The orbital period 



Previous results (Niemela 1991) already have shown the variability of radial 
velocities of the emission lines in the spectrum of WR 98. We searched for 
periodicities of the radial velocity variations in the strongest emission lines in 
our spectra, namely He A 4686A and Cm A 4648A using algorithms published 
by Marraco & Muzzio (1980) and Cincotta et al. (1995). The most probable 
period obtained by both codes is 47.8 days. 

Using this period as initial value, we then calculated orbital solutions for 
the radial velocity variations of both emission lines. We found that the so- 
lutions tend to circular orbits, with the most probable value of the period 
as 47 .825days ± 0.005. However, other periods close to this value can not be 
discarded. 
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4000 4400 4800 

Wavelength [A] 



Fig. 2. Continuum normalized spectra of WR123 (WN8o), WR98, and WR55 
(WN7o). Note the similarity of the WN spectrum of WR 98 with that of WR 55. 

3.3 The radial velocity orbit 

Orbital elements for each emission line and the H7 and H/3 absorptions were 
determined with an improved version of the program originally published by 
Bertiau & Grobben (1968). With the present data, the orbits of the four 
emission lines, Niv, Nv, Cm, and Hell have negligible eccentricity, thus we 
have fitted circular orbits for all our radial velocities. Circular orbital elements 
for the individual lines are listed in Table 3.3, where Vq refers to the center-of- 
mass velocity, K, to the semi-amplitude of the radial velocity variations, and 
T is the time when the WR star is in the front of the system. 

In our orbital fits the three best defined emission lines gave equal T values 
within errors, therefore we adopt as ephemeris for the WR 98 binary system: 

T = 2,445,676.4+47.825E 

Niv, Nv, Cm, and Hell emission lines move in phase, indicating that they 
are formed in the same stellar envelope. This is also the case in the two other 
known WN/C binaries, namely WR 145 (MR 111) and WR 153 (GP Cep) 
(Massey & Grove 1989). No detectable phase delays among emission lines are 
present within the errors of the orbital fits. 

Semi-amplitudes of the orbital motion of the Hell and Cm emission lines 
appear lower than those of ionized Nitrogen emission lines. This effect is also 
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4830 4860 4890 

Wavelength (A) 



Fig. 3. Continuum rectified spectra of the Hen 4859 A emission with the superim- 
posed H/3 absorption observed during four different phases of the binary system. 
Note the antiphased movement of absorption and emission. The spectra are inten- 
sity-shifted for a better presentation. 

observed in other WR binaries, e.g. WR 29 (Niemela & Gamen 2000), and 
may arise if the Hell and Cm lines are partly formed in the interaction region 
of the binary components. 

We measured H7 and H/3 absorptions in those spectra of WR 98 were this 
was possible. Radial velocities of these hydrogen lines phased with the binary 
period move anti-phased with the emission lines, thus indicating that they 
belong to an O type companion of the WR component in the binary system. 
Fig. 3 depicts the behaviour of the H(3 absorption line upon the Hell 4859 A 
emission in four different binary phases, illustrating the antiphased movement 
of the absorption and emission lines. 

The circular orbital elements for the hydrogen absorption lines are included 
in Table 3.3, and depicted in Fig 4. 

The orbital semi-amplitudes of the radial velocity variations of the absorp- 
tion lines and of Niv, Nv, Hell and Cm emission lines indicate mass-ratios 
(•M-wr/ ' -M-o — <?) between 1 and 1.8. If a monotonic outward decreasing tem- 
perature gradient exists in the expanding WR envelope, we expect that the 
highest ionization emission, namely Nv, represents better the orbital motion 
of the WN/C component, thus a value of q ~ 1 seems more plausible. How- 
ever, taking into account the rather high uncertainties in the radial velocity 
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Table 2 

Circular Orbital Elements of WR 98 



WN/C OB 
Parameter N IV em. N v em. C in em. He n em. absorptions 



P [days] 47.825 ± 0.005 



Vblkms- 1 ] 


-41 ±4 


-15 ±3 


-47 ±2 


4±2 


6 ±15 


Kfkms- 1 ] 


106 ±6 


109 ±5 


72±3 


65 ±3 


112 ±16 


T [HJD]W 


6.4 ±0.5 


5.5±0.4 


6.5±0.3 


6.3 ±0.3 




a sin i [R Q ] 


100 ±5 


103 ±5 


68 ±3 


61 ±3 


106 ±15 


Mwn sin 3 i [Mq] 


27 ±10 


28 ±10 


19 ±8 


18 ±7 




Mob sin 3 i [Mq] 


25 ±7 


27±7 


12 ±3 


10 ±3 




Q 


1.06 


1.03 


1.56 


1.72 





* HJD 2,445, 670±: Time of conjunction, the WN star in front of the system. 



values of the absorption lines, these mass-ratios should not be overinterpreted. 
Values of minimum masses and q for both components are tabulated for each 
emission line in Table 3.3. 

We also estimated the radii of the critical Roche lobes using the expression 
given by Paczynski (1971), which resulted r RL sin i = 79.5i? Q for both com- 
ponents. Each component of the WR 98 binary system seems to be well inside 
its critical Roche radius. 

3.4 The radial velocity analysis of P-Cygni absorption lines 

We analysed the radial velocities of the three strongest P-Cygni absorption 
lines in our spectra of WR 98, namely Hel 3888A, 447lA, and 5015A. As the 
Hel 3888A absorption presents a rather asymmetric profile, we measured the 
barycenter of this line. Because Hel 50 15 A was not observed in the wavelengh 
range of the photografic spectra, the radial velocities of this line were only 
determined in the digital spectra. The three P-Cygni absorption lines follow 
the orbital motion of the WN/C component of the binary, as illustrated in 
Fig 5. 

Circular orbits were fitted to the radial velocity variations of the P-Cyg ab- 
sorptions. We obtained similar systemic velocities for the three absorption 
lines: v ~ —1130 ± \hkms~ l . This systemic velocity could be considered as 
a lower limit of the terminal velocity of the WR stellar wind, and in fact the 
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0.5 1 1.5 

Orbital phase 



Fig. 4. Radial Velocities of Hen, Cm, Niv, and Nv emission lines, and H7 (circles) 
and H/3 (triangles) absorption lines phased with the ephemeris 2,445, 676. 4+47.825E. 
Open symbols represent the photographic data, and solid symbols, the digital CCD 
data. Curves are the orbital solutions from Table 3.3. 
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Fig. 5. Variations of the radial velocities of two Hel P-cygni absorption lines in the 
spectrum of WR 98, symbolized and phased with the same ephemeris as emission 
lines in Fig. 4. Dashed lines represent circular orbital solutions. 

value agrees with the determination of terminal wind velocity of WR 98 by 
Eenens & Williams (1994). 

Semi-amplitudes of the radial velocity variations of the P-Cyg absorptions 
gave different values, namely 65 and 104 kms^ 1 for Hel AA 3888, and 4471 
and 5015A respectively. The lower semi-amplitude of the orbital motion of 
Hel 3888 A could be indicating that this line has another component, possi- 
bly originating in a common expanding envelope surrounding the binary, and 
which is not posible to deblend in our spectra. Similar behavior of Hel A 3888A 
low-energy metastable absorption line is observed in the spectra of other WR 
binary stars, e.g. HD 214419 (= WR 155) (See Leung et al. 1983). 



4 Conclusions 



Our radial velocity analysis of lines in the spectrum of WR 98 based on a 
long term spectroscopic database confirms that this star is a binary system 
with a period of 47.825 days. Nitrogen and Carbon emission lines in the spec- 
trum move in phase, which indicates that they are formed in the same stellar 
envelope. 

We have detected faint absorption lines upon the WR emissions in our spectra. 
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These absorptions belong to an O-type component of the binary, with orbital 
motion anti-phased with the WR emissions. Thus, WR 98 is one more member 
of the very limited sample of WR+OB double-lined binaries, the first one with 
a WR component which appears to be at an intermediate evolutionary phase 
between the WN and WC stages. 

The minimun masses indicated by our radial velocity orbit are quite high. 
Adopting the radial velocity orbit of the Nv emission as representative of the 
motion of the WR component of the binary, we obtain values of minimum 
masses of 28 and 27 M for the WR and OB components, respectively. There- 
fore, the component which at present shows the WR type spectrum, must have 
been very massive during its main sequence stage. This is in agreement with 
the estimate of 80 M for the cluster turn-off of Tr 27 (Massey et al. 2001), 
of which WR 98 is a probable member. 

A comparison of the WN spectrum of WR 98 with other stars classified as 
WN8o (WR 123) and WN7o (WR 55), shows a remarkable resemblance of 
WR 98 with the WN7o spectrum, as seen in the Fig. 3. In our data, the 
classification of WR 98 as WN7/C appears more likely. 

P-Cygni absorption lines of Hel in our spectra of WR 98 also follow the orbital 
motion of the WR component of the binary, thus these lines are formed in the 
expanding atmosphere of the WR star. We determined a systemic velocity 
of the P-Cyg absorptions of ~ — 1130A;m for WR 98, which is in good 
agreement with the value of the terminal wind velocity published by Eenens 
& Williams (1994). 
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